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Abstract. We investigate the need and prospects for measuring dark matter properties at particle collider experiments. We 
discuss the connections between the inferred properties of particle dark matter and the physics that is expected to be uncovered 
by the Large Hadron Collider (LHC) and the International Linear Collider (ILC) and motivate the necessity of measuring 
detailed dark matter properties at a collider. We then investigate a model-independent signature of dark matter at a collider 
and discuss its observability. We next examine the prospects for making precise measurements of dark matter properties using 
two example points in minimal supergravity (mSUGRA) parameter space. One of the primary difficulties encountered in such 
measurements is lack of constraint on the masses of unobservable heavy states. We discuss a new method for experimentally 
deriving estimates for such heavy masses and then conclude. 
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INTRODUCTION 

The existence of non-baryonic, non-luminous matter ap- 
pears now to be firmly established. The amount of such 
'dark' matter has been measured at the 10% level. The 
current best limit is given in terms of the fractional 
energy density of the universe that is contained in dark 
matter, and h, the reduced Hubble parameter 1 1 ]: 

0.094 < Q. x h 2 < 0.129 (at 2a). (1) 

The nature of dark matter is one of the greatest un- 
solved mysteries of science. Usually, one hypothesizes 
an extra stable electrically neutral particle that can serve 
as dark matter. Most particle physics models that at- 
tempt to explain the weak scale and electroweak sym- 
metry breaking contain such additional particles among 
their non-Standard Model particle spectra. Since one of 
these models might be discovered at the LHC or ILC, it 
is interesting to ask what, if anything, such colliders can 
tell us about dark matter. This talk focuses specifically 
on answers to this question. 

CONNECTIONS BETWEEN DARK 
MATTER AND COLLIDERS 

Precious little is known about dark matter, except for its 
gross density, given in Eqn. ^ If one assumes that dark 
matter comes from a thermal relic, this gross density is, 
in large part, determined by the dark matter annihilation 



1 Plenary talk given at PASCOS05, Gyeongju, Republic of Korea, June 
2005. 



cross section: o(%% — > SMSM). In fact, the present- 
day dark matter abundance is roughly inversely propor- 
tional to the thermally averaged annihilation cross sec- 
tion times velocity 2 : i2 z /i 2 «1/ (<7v). Figure^shows the 
constraint on the annihilation cross section as a function 
of dark matter mass that results from Eqn. ^ Remark- 
ably, this constraint is effectively independent of dark 
matter mass 3 , and points to cross sections expected from 
weak-scale interactions (around 0.8 pb for s-annihilators 
and 6 pb for /5-annihilators). It points to the possibility 
that the solution to the dark matter puzzle is connected to 
the explanation for the weak scale, leading to the idea of 
Weakly Interacting Massive Particle (WIMP) dark mat- 
ter. Such WIMPs exist not only in supersymmetric the- 
ories, but also in theories involving extra dimensions |2] 
and 'little Higgs' theories |3] as well. The LHC and the 
ILC are specifically designed to probe the origin of the 
weak scale, so dark matter and future collider physics 
appear to be intimately connected. 

Many experiments exist that hope to discover dark 
matter at many distances from earth. However, they, 
by themselves, are not able to solve the dark matter 
puzzle. Direct detection experiments hope to measure 
the mass of dark matter particles and their scattering 
cross section off atomic nuclei, o{%N^%N). Since 
any measured rate of scattering depends on the num- 



2 Usually the cross section times velocity can be conveniently ex- 
panded in powers of relative dark matter particle velocity: av = 

Y.J Gem v^ 2J K Commonly only the lowest order non-negligible power 
of v dominates. For 7 = 0, such dark matter particles are called s- 
annihilators, and for J = l, they are called p-annihilators; powers of 
J larger than l are rarely needed. 

3 This effect is due to the changing number of degrees of freedom at 
the time of freeze-out as the dark matter mass is changed. 
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FIGURE 1. Values of the quantity a an allowed at 2o" level 
as a function of the dark matter mass. The lower (upper) band 
is for models where s-wave (p-wave) annihilation dominates. 
Taken from Ref. jj. 



ber of dark matter particles present, astrophysical uncer- 
tainties in halo models affect the accuracy of measure- 
ments. Furthermore, to claim that the dark matter puzzle 
has been solved, we must determine that a given par- 
ticle (or set of particles) can provide the correct (<7v) 
shown in Fig. ^ It is generically not possible to relate 
o(%N %N) to the required a {%% — > SMSM). Indi- 
rect detection experiments attempt to discover dark mat- 
ter through its present-day annihilation in regions of high 
density such as the galactic center. Again, significant as- 
trophysical uncertainties exist. Also, these experiments 
hope to measure either a single annihilation final state, 
such as <7 — » v v )> or a subset of final states, such as 
those with charged particles in the final state. While these 
final states do contribute to a (%% — > SMSM), again, 
there is not a model-independent relation between the 
two. As a result, it is necessary to make detailed mea- 
surements of the particles and couplings involved in dark 
matter annihilation. These detailed measurements can 
only be done at particle collider experiments, such as the 
LHC and ILC. 



MODEL-INDEPENDENT COLLIDER 
SEARCH FOR DARK MATTER 

Since the WIMP hypothesis and the known abundance 
of dark matter give specific values for the dark matter 
annihilation cross section, one might hope this cross 
section can be turned into a rate for a measurable process 
at a particle collider. Such a model-independent collider 



signature does indeed exist |4] 4 and has implications for 
indirect detection experiments 1 5 ] . 

To arrive at the model-independent signature, we start 
with cosmological data, namely the allowed values for 
a un from Fig. [2 We introduce the parameter K e : 

K e = a(xx^e + e-)/a{xx^SMSM) (2) 

This equation relates the dark matter annihilation 
cross section to the cross section involving e + e~ in 
the final state. Then we use crossing symmetries to re- 
late (7 (XX ~* e + e~) to c (<? + e~ — > %%)■ Finally, we use 
collinear factorization to relate o(e + e~ — > %%) to the 
process o(e + e~ — ► XX7)- We thus have gone from cos- 
mological data, in the form of rj„„, to a prediction for the 
rate e + e~ — > XX 7'- 
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Here x = 2Ey/^fs, 9 is the angle between the photon 
and the incoming electron, S x is the spin of the WIMP, 
and Jq is the dominant value of J in the velocity expan- 
sion for (Tv, defined earlier. 

The accuracy of our method is illustrated in Figure |2 
The left panel in Figure |2 shows how the formula in 
Eqn.|5] (green line, labelled as 'Eq. (9)') compares with 
the exact calculation in a supersymmetric theory (red 
line, labelled 'MSSM') with a WIMP mass of 225 GeV. 
It is clear that the agreement is quite good, especially 
for the harder photons, which are most useful in our 
approach. The right panel shows the expected reach in 
K e for a 500 GeV linear e + e~ collider as a function of 
the WIMP mass. We have included background from 
e + e~ — > wy. We show the 3 <T contour (black solid 
line) including statistical errors and also indicate the 
values one might expect from supersymmetric models 
(the region below the red dashed line labelled 'SUSY'). It 
is clear that this signal, while challenging, could provide 
direct evidence of WIMP production at colliders. 



TESTING SUPERSYMMETRIC 
COSMOLOGY AT THE ILC 

The above model-independent rate is rarely the domi- 
nant collider signature of new physics within a given 



4 Due to space restrictions, here we are skipping many intermediate 
steps. For further detail, we refer the reader to the original paper, 
Ref. U. 
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FIGURE 2. Left panel: Comparison between the photon spectra from the process e + e~ — * 2Xi + y in the explicit supersymmetric 
models defined in Ref.|4] (red/dark-gray) and the spectra predicted by Eqn.|3|for a p-annihilator of the corresponding mass and 
Ke (green/light-gray). Taken from Ref. |4]. Right panel: The reach of a 500 GeV unpolarized electron-positron collider with an 
integrated luminosity of 500 fb _I for the discovery of p-annihilator WIMPs, as a function of the WIMP mass M% and the e + e~ 
annihilation fraction K e . The 3 0" (black) contour is shown, along with an indication of values one might expect from supersymmetric 
models (red dashed line, labelled 'SUSY'). Only statistical uncertainty is included. 



model. It therefore makes sense to look at more model- 
dependent processes; we choose the minimal supersym- 
metric extension of the Standard Model (MSSM). The 
WIMP is taken to be the lightest supersymmetric parti- 
cle (LSP), usually the lightest neutralino, %\ - If tne LSP 
is the dark matter particle, the generic collider signatures 
of supersymmetry all involve missing energy due to the 
two stable Xi s escaping the detector. In order to prove 
that the missing energy particle is indeed a viable WIMP 
dark matter candidate, one needs to calculate its expected 
present relic abundance. One thus needs to measure all 
parameters which enter the calculation of (<7v). In the 
most general MSSM there are more than 100 input pa- 
rameters at the weak scale, but fortunately, a lot of them 
are either tightly constrained or not very relevant for the 
dark matter calculation. Nevertheless, there are still sev- 
eral relevant parameters left; they must be determined 
from collider data. Of course, the relevance of any one 
parameter depends sensitively on the parameter space 
point. We now look at two points: first, a point from the 
'bulk region,' and then a point from the 'focus point re- 
gion,' [6] both from mSUGRA, a convenient subset of 
the larger MSSM 0110]. 

Our first point, point B' from the bulk region, is de- 
fined by the following input parameters 5 : mo — 57 GeV, 
m 1/2 = 250 GeV, A Q = 0, tan/3 = 10 and sign(ju) = +1. 
Fig. [3] shows that the particle spectrum at this point 
is quite light. The two lightest neutralinos, the light- 



Here we have used Isajet 7.69 0^1 and DarkSUSY [TJ except for 
calculations where coannihilations are important, in which case we 
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FIGURE 3. Mass spectrum of point B' from the bulk region. 
Taken from Ref. 191. 



est chargino and all of the sleptons have masses below 
200 GeV. All of the squarks are lighter than 600 GeV. 
The heaviest particle, the gluino, only weighs 611 GeV. 
Therefore, one would expect colliders to have significant 
discovery and measurement capabilities. 

We first estimate the sensitivity of £l x h 2 to the various 
SUSY parameters, and then we determine the precision 
with which they can be measured at colliders. In Fig.0] 
we show the sensitivity of the dark matter relic density 
to 6 relevant MSSM parameters. In each panel, the green 
region denotes the 2a WMAP limits on £l x h 2 and the 
red line shows the variation of the relic density as a func- 
tion of the corresponding parameter. The vertical (blue- 
shaded) bands denote parameter regions currently ruled 
out by experiment. The blue dot in each panel denotes the 



nominal value for the corresponding parameter at point 
B'. 

The importance of the different SUSY parameters for 
the determination of the neutralino relic density can be 
judged from the slope of the lines in Fig. 0] If £l%h 2 is 
insensitive to a given parameter, the corresponding line 
will be flat. Conversely, if the relic density is particularly 
sensitiive to some SUSY parameter, the slope of the cor- 
responding variation curve will be very steep. Interest- 
ingly enough, Fig. |4] also shows that even for one given 
parameter, there are regions where this parameter can be 
important (for example Mj }r at low masses, where coan- 
nihilation processes with bottom squarks become impor- 
tant), as well as regions where this parameter has very 
little impact on the relic density (e.g. Mg at its nom- 
inal value). Therefore, in estimating the uncertainty on 
£l z /z 2 , collider data on otherwise "irrelevant" parameters 
can be very important. In the absence of any informa- 
tion about the value of a given parameter, one should let 
it vary within the whole allowed range, which may en- 
compass values of the parameter for which it becomes 
relevant. 

The behavior of some of the lines in Fig. |4] can be 
understood as follows. At point B' the lightest neutralino 
%i is mostly Bino, hence one would expect that its relic 
density will be sensitive to the Bino mass parameter M\ . 
Indeed, this is confirmed by Fig.|4ja). For small M\, we 
observe enhanced sensitivity near the Z and Higgs pole 
regions (2M% ~ Mz and 2Mj> ~ M/,). For large values of 
Mj we see a significant variation again, this time because 
the neutralino LSP becomes more and more degenerate 
with the sleptons, and its relic density is depleted due 
to coannihilation processes. The analysis of Figs. |4jb) 
and |4]c) is very similar: the sfermions are irrelevant, 
if they are heavy, but may become very important if 
they are sufficiently light to induce coannihilations. For 
explanations of the behavior of the other lines, please 
refer to Ref. £J. 

Having determined the correlations between the weak- 
scale SUSY parameters and the relic abundance of neu- 
tralinos, it is now straightforward to estimate the uncer- 
tainty in Q. x h 2 after measurements at different colliders. 
The result is shown in Fig.|5J where the outer red (inner 
blue) rectangle indicates the expected uncertainty at the 
LHC (ILC) with respect to the mass M x and relic density 
£l z /z 2 of the lightest neutralino. The yellow dot denotes 
the actual values of M x and £l x h 2 for point B' and the 
horizontal green shaded region is the current measure- 
ment, Eqn.n 

In arriving at this result, we made the following as- 
sumptions about measurements at the LHC. We assume 
measurements of the jjj^, Xo, an d Xi masses at the level 
of 10%. The precision on left-handed squark masses 
should be no better than the precision on gaugino masses, 



but we have assumed 10% again. We have conserva- 
tively assumed that no slepton or right-handed squark in- 
formation will be available. Finally, in terms of Higgs 
bosons, we expect a detection only of the lightest (Stan- 
dard Model-like) Higgs boson and the absence of a heavy 
Higgs boson signal will simply place the bound > 
200 GeV. 

Our assumptions about the corresponding precision at 
a 500 GeV ILC were the following. Since superpartners 
need to be pair-produced, we take all sparticles lighter 
than 250 GeV to be observable, and their masses can be 
measured to within 2%. This includes the same chargino- 
neutralino states as in the case of LHC, plus all sleptons. 

From Fig.|5]we see that with the assumptions above, 
the LHC (scheduled to turn on in 2007) is not competi- 
tive with the current state of the art determination of the 
relic density from cosmology. Nevertheless, it will bound 
the relic density from above and below, and may provide 
the first hint on whether the dark matter candidate being 
discovered at colliders is indeed the dark matter of cos- 
mology. The ILC will fare much better, and will achieve 
a precision rivalling that of the cosmological determina- 
tions. 

Turning to point LCC2, we show in Fig. [5] the spec- 
trum of particle masses. We also show the analogous 
variation of the relic density as a function of 5 rele- 
vant parameters. Here the squark and slepton masses are 
heavy and have little impact on the actual £l%h 2 . For this 
point, the lightest neutralino is mostly Bino, but with a 
non-negligible higgsino component. 

Fig-EJa) shows the dependence of Q. x h 2 on the Bino 
mass parameter Mi . As expected, lowering Mi increases 
the Bino component of the LSP, thus suppressing (7 an 
and increasing ii^/i 2 . Fig. |6jb) exhibits complementary 
behavior: the M2 parameter controls the wino fraction of 
the LSP, and small values of M2 lead to wino-like dark 
matter, which has a large annihilation rate and therefore 
smaller relic abundance . Again, see Ref. 1 9] for further 
details. 

These results will be combined with the outcome of 
a comprehensive simulation study, including detailed de- 
tector simulation, on the expected experimental precision 
at a 500 GeV ILC for point LCC2. The final product will 
be the analogue of Fig.|5]for point LCC2. For further de- 
tails on the current status of this analysis, see 0. 



6 Neutralino dark matter with significant wino content has been studied 
in the contexts of scenarios with non-universal gaugino masses 1 1 311 _ 
string-derived supergravity 1 14] and anomaly-mediated supersymmetry 
breaking 11 II . 





100 200 300 400 500 

M T (GeV) 



100 200 300 500 700 1000 2000 3000 

'■A ( GeV > 




tan0 



FIGURE 4. Effect on relic density of varying relevant SUSY mass parameters for point B'. The horizontal (green-shaded) region 
denotes the 2a WMAP limits on the dark matter relic density. The red line shows the variation of the relic density as a function of 
the corresponding SUSY parameter. The vertical (blue-shaded) bands denote parameter regions currently ruled out by experiment. 
The blue dot in each plot denotes the nominal value for the corresponding parameter at point B'. Taken from Ref. |9]. 
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FIGURE 5. Accuracy of WMAP (horizontal green shaded 
region), LHC (outer red rectangle) and ILC (inner blue rectan- 
gle) in determining M%, the mass of the lightest neutralino, and 
its relic density Q.%h 2 . The yellow dot denotes the actual values 
of M x and Q x h 2 for point B\ Taken from Ref. |jj and Ref. Q|. 

MEASURING SLEPTONS AT THE LHC 

As discussed in the previous section, slepton masses are 
among the key parameters in determining whether %\ is a 
good dark matter candidate 1 16]. Without a collider mea- 
surement of the slepton mass, there may be a significant 
uncertainty in the relic abundance calculation. This un- 



certainty results because the slepton mass should then be 
allowed to vary within the whole experimentally allowed 
range. 

The importance of slepton discovery is two-fold. First, 
supersymmetry predicts a superpartner for every stan- 
dard model particle. Therefore, the discovery of the su- 
perpartners of the leptons is an important step in ver- 
ifying supersymmetry. Second, knowledge of slepton 
masses is always important for an accurate determina- 
tion of the relic abundance of %y. 

Here we show that the LHC will indeed have sensi- 
tivity to mSUGRA slepton masses, even in the case of 
heavy sleptons fl7ll . We show that the difference between 
real and virtual sleptons can be clearly seen. We also 
show that in case of virtual sleptons, one can frequently 
limit the allowed range of their masses, which is equiva- 
lent to a rough indirect slepton mass measurement. 

Slepton studies at the LHC are challenging. Direct 
production of sleptons suffers from large backgrounds, 
mostly due to and tt production 1 18]; the meth- 

ods for slepton mass determination used at the linear col- 
lider are not applicable here. Fortunately, sleptons are 
produced in sizable quantities at the LHC through cas- 
cade decays. These events can be easily triggered on 
and separated from the standard model backgrounds. A 
common hierarchy in supersymmetric models is \M\ \ < 
IM2I < |ju|. In that case, sleptons affect the decay — > 
£^£*Xi- The resulting dilepton distribution, in principle, 
contains information about the slepton mass mi. This sit- 




uation is complicated by the fact that %^ can also decay 
through a real or virtual Z: — ► zMjfj 5 — > ^ ± ^ T ^?- 

What is the observable in these events that is sensitive 
to the slepton mass? In this analysis we will consider the 
dilepton invariant mass distribution, mu- The endpoint 
of mu contains information about the masses of the real 
particles involved in the decay But there is more 
information contained in the mu distribution than just in 
the endpoint ifTvIl - One would expect the shapes of the 
Z and £ mediated distributions to be different. Further- 
more, the shape of the total decay distribution (including 
both Z and t contributions) should change as a function 
of the slepton mass. Indeed this is the case, as illustrated 
in Figure [7] Here we show the dilepton invariant mass 
distribution resulting from the interference of the Z and 
eg -mediated diagrams. The kinematic endpoint is kept 
fixed; this illustrates that endpoint analyses are largely 
insensitive to the slepton mass. In all four cases, the slep- 
ton is virtual, but we see a clear difference in the shape 
of the distribution. In the case of two body decay through 
a real slepton, the mu distribution will be triangular, as 
required by phase space. It is also clear that the mu dis- 
tribution changes significantly as a function of slepton 
mass 

How does one use this fact to statistically differenti- 
ate between distributions coming from sleptons of differ- 
ent mass? The Kolmogorov-Smirnov (K-S) test provides 
a statistical procedure for differentiating between such 
distributions [20]. The K-S test calculates the maximum 
cumulative deviation between two unit-normalized dis- 
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FIGURE 7. Comparison of M e + e - distributions for different 
selectron masses. We only consider the diagrams mediated by 
Z and en. All parameters are held fixed except for the en mass. 
The (green, blue, red, black) line is for a (300, 500, 1000 GeV, 
and oo) mass selectron, respectively. The neutralino masses 
are kept constant, and their difference is 88 GeV. Taken from 
Ref. El. 

tributions. It then translates this number into a confidence 
level with which it can be excluded that the two distribu- 
tions come from the same underlying distribution. So, the 
K-S test cannot tell that two distributions came from the 
same underlying distribution, but it can tell that two dis- 
tributions did not come from the same underlying distri- 
bution. We use this test to explicitly investigate the ability 
of the LHC to determine slepton masses in mSUGRA. 




M„(GeV) 



FIGURE 8. Slepton mass determination in mSUGRA with 
A(> = and tan/3 = 10. The plot shows the effect on the 
mSUGRA parameter space of fixing the dilepton kinematic 
endpoint of the — > e + e~Xi decay to be mumax — 59 GeV. 
It also shows the ability of the K-S test to bound the mass of 
the lightest selectron for light virtual selectrons. Details are 
explained in the text and Ref. fl7ll . 



We now demonstrate a slepton mass measurement at 
the LHC in mSUGRA parameter space. We select 1000 
events that contain the decay of %2 ~* e+e X\ - We imag- 
ine that the LHC experiments have observed the dilepton 
mass distribution and have measured a kinematic end- 
point at 59 GeV. The points with this same kinematic 
endpoint are the solid lines (blue for > 0, black for 
H < 0) in Figure [8] The dashed lines in the upper left- 
hand corner indicate where m^o = m?, . The dashed lines 
running through the middle of the plot indicate where 
mg R — m^o- This is where the slepton-mediated neu- 
tralino decays changes from being three body {eg is vir- 
tual to the right of these lines) to two body (Sr is real to 
the left of these lines). 

In Figure [8] we also display the power of the K-S test 
in identifying the mass range of the slepton. We have 
taken one point (denoted by the black star) as result- 
ing from a possible experimental measurement involving 
1000 signal events. Then we have compared this 'exper- 
imental data' with template distributions (denoted by the 
red and yellow dots). The red dots indicate points that the 
K-S test can identify as not coming from the same distri- 
bution as our 'experimental data' at the 95% confidence 
level. The yellow dots denote points that could not be 
excluded at the 95% confidence level. We can rule out at 
the 95% confidence level almost all of the example tem- 
plate points except for the ones with the same value of ji, 
but with almost identical values of Mq. Thus, the mass of 
the slepton can be determined with some accuracy even 
though it is not produced as a real particle. Analysis of 
experimental points with other values of the slepton mass 
can be found in Refs. 1 171 12111 . 



CONCLUSION 

In this talk, we have introduced techniques and concepts 
for measuring dark matter at colliders. We have moti- 
vated why this should be possible and also why it is ab- 
solutely necessary to validate any claim of a solution to 
the dark matter puzzle. We have discussed a challeng- 
ing model-independent signature of dark matter. We have 
also investigated the prospects for making precise mea- 
surements of cosmologically relevant parameters in a 
model-dependent approach. Finally, we discussed a new 
technique with implications for dark matter determina- 
tions at a collider that allows slepton masses to be mea- 
sured at the LHC, even for heavy virtual sleptons. 
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